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We fit the HiRes uhra-high energy cosmic ray (UHECR) spectrum measurements with broken power laws in 
order to identify features. These fits find the previously observed feature known as the Ankle at 10^*'^ eV, 
as well as evidence for a suppression at higher energies, above 10^^'^ eV. We use the integral spectrum and 
the Ei^^ test to identify this high energy suppression with the GZK suppression. Finally, we use a model 
of uniformly distributed extragalactic proton sources together with a phenomenological model of the galactic 
cosmic ray spectrum to compare the HiRes spectra to what should be expected from the GZK suppression, and 
to measure how the extragalactic sources must evolve and what the input spectral slope must be to fit the HiRes 
data. Fits using updated spectra will be presented in Pune. 



1. Broken Power Law Fits and 

The HiRes Collaboration has recently released two measurements of the UHECR flux using monocular ob- 
servations from its two sites] 1 1. Broken power law (BPL) fits made to these measurements can be used to 
identify features in the spectrum and to estimate their statistical significance in the data. All the fits presented 
in this paper are performed using the normalized binned maximum likelihood method] 2 1. This method re- 
quires comparing the numbers of events expected in a given model to the number actually observed. The 
number of expected events is obtained from the predicted flux divided by the same exposure used to calculate 
the observed flux. The numbers of events observed by HiRes are shown in Figure ^ The binned maximum 
likelihood method also allows one to use bins in which there were no observed events, but in which events were 
expected. The result of the fits are expressed in terms of a quality-of-fit parameter which approaches a true 
in the limit of large numbers of events. The BPL fits are only made to bins above 10^^^ eV. While HiRes-II 
does have three bins below this energy, these bins have poor statistics. We don't fit these bins in order to avoid 
biases due to an expected change in the spectral slope, the Second Knee, at about this energy. Parameters from 
fits to the HiRes monocular spectra to a BPL, J{E) = CE with zero, one and two floating break points are 
shown in Tabled The result of the fit in the two floating break point case is shown in Figure|2] 

The simple power law fit is clearly not a good fit. Adding one floating break point gives a much better fit, 
and the break point finds the feature known as the Ankle with a very high degree of statistical significance. 
Adding a second floating break point improves the fit further The break point is found to be at approximately 
the energy expected of the GZK suppression] 3 1. 

The statistical significance of the second break in the spectrum can be estimated by looking at the reduction in 
the achieved by adding the break point, or by comparing the number of expected events above the second 
break point to what would be expected if the spectrum continued unabated above the second break point. The 



Table 1. Parameters found in broken power law fits to the HiRes monocular spectra. 
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Figure 1. The numbers of events in each bin of 
the HiRes monocular spectrum measurements HI). The 
HiRes-I measurement includes two empty bins cen- 
tered at 20.3 and 20.5. The HiRes-II measurement in- 
cludes two empty bins centered at 20.1 and 20.3. 




Figure 2. The HiRes spectra fit to a BPL with two 
break points. The parameters of the fit are given in Ta- 
bleQ The red line is used to calculate the significance 
of the second break point and to calculate the expected 
integral spectrum in the Ei^^ calculation. 



reduction of 16 in the for two additional degrees of freedom corresponds to just under 4cr significance in a 
gaussian fit. Using the red line in Figure|2] one would expect to see 28.0 events, where 1 1 events were actually 
observed. The Poisson probability of observing 11 events or fewer when expecting 28.0 is 2.4 x 10"'*. For 
comparison, the area in one tail of a gaussian distribution outside of 4a (3a) is 3.2 x 10~^ (1.4 x 10"'^). So 
the significance of the high energy suppression is between 3cr and 4a. 

Berezinsky et a/.0 have suggested measuring the energy of a break in the UHECR spectrum finding the 
energy, Ei^^, at which the observed integral spectrum is half of what one would expect with no break. The 
integral spectrum measured by HiRes, along with the expected integral spectrum using the red line in Figure|2l 
is shown in Figure |3] The ratio of the observed to the expected integral spectra is shown in Figure |3 By 
interpolating between the HiRes-1 points in Figure|3 we find an experimental value of logj^Q Ei^ ~ 19.77^q Jg 
{E in eV). Berezinsky et al.\4j have determined logio ^i/2 ~ 19-'''2 as what is to be expected for the GZK 
suppression for 2.1 < 7 < 2.7. 



2. Uniform Source Model Fits 

Berezinsky et al. [4j also calculate the energy loss rate for UHE protons traveling through the cosmic microwave 
background radiation, assuming continuous energy loss through electron pair production, pion production and 
universal expansion. This can be used to predict the observed spectrum of extragalactic protons for a given 
input spectrum. 

Since protons with energies above the pion production threshold lose energy in highly inelastic interactions, we 
have used the Monte Carlo method of DeMarco et al.^ to model this part of the propagation of extragalactic 
protons. Protons from a shell at a given redshift are propagated from generation to observation (at z = 0). 
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Figure 3. The integral spectra derived from thie dif- 
ferential spectra shown in Figure |2| The blue points 
represent the expected integral spectrum from the red 
line in that figure. 



Figure 4. The observed-to-expected ratio for integral 
spectra. The value of Ei^^ is obtained by a simple in- 
terpolation between the HiRes-I points as indicated by 
the blue lines. 



Each input energy corresponds to a distribution of observed energies after such a propagation. This Eii -Eo^t 
distribution convolved with the input spectrum gives the observed spectrum for protons from sources at a given 
redshift. If the distribution of sources is assumed to be uniform at any given redshift, but to evolve with redshift 
as (1 + z)™, one can combine the spectra from different shells as shown in for a coarse, logarithmic series of 
shells in Figure|5] 

Because the HiRes-II spectrum extends to fairly low energies, one should take into account an additional 
galactic component when fitting to this USM model. We made the simple phenomenological assumption that 
the extragalactic and galactic components of the spectrum are expressed, respectively, in the light (protons) 
and heavy (iron) components of a composition measurement. We use fits to the HiRes Prototype/MIA|6| and 
HiRes Stereo|7 1 composition measurements with respect to QGSJet proton and iron expectations to determine 
the relative sizes of the extragalactic and galactic components. Then, for a given set of parameters for the 
extragalactic UHECR spectrum, we can add the appropriate galactic cosmic ray flux. 

We varied the input spectral slope, 7, and evolution parameter, m, to find the best fit of the HiRes monocular 
data to this USM-plus-Galactic model. The best fit spectrum is shown in Figure |6l The final result for the 
extragalactic USM model is 7 = 2.38±0.035(stat)±0.03(syst), m ^ 2.55 ± 0.25(stat) ±0.30(syst), where 
the systematinc uncertainties come from varying the extragalactic/galactic ratio within the limit allowed by the 
composition measurements. 

The fit works best in the region on either side of the Ankle, with the fall at lower energies, into the ankle, 
primarily determining m, and the rise at higher energies, out of the Ankle, primarily determining 7. The fit 
does not work as well in the region just below the GZK suppression. This could be due to some sources having 
a maximum energy below the GZK threshold. The is also little sign of a Second Knee at around 10^^'^ eV, 
though the HiRes-II data has little statistical power in this region. 
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Figure 5. Spectrum of extragalactic protons from 
USM with 7 — 2.4 and m — 2.5. The spectra indi- 
vidual shells sum to the red line shown. A finer series 
of shells sum to the black line. 



Figure 6. Best fit of the USM-plus-Galactic model 
to the HiRes monocular spectra. The red line shown 
the extragalactic component, the green line shows the 
galactic component and the black line shows the sum. 



3. Conclusion 



The HiRes detector has observed the Ankle and has evidence for a suppression at higher energies above 10 
eV. The energy for this high energy suppression agrees with what is expected from the GZK suppression 
according to the Ei^^ test. The observed spectra are well fit by a USM-plus-Galactic model, which finds 
an input spectral slope for extragalactic protons of 7 = 2.38 ± 0.035(stat) ± 0.03(syst), and an evolution 
parameter m = 2.55 ± 0.25(stat) ± 0.30(syst). 
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